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ABSTRACT 

Let  Mjj  denote  the  bivariate  box-spline  corresponding  to  the  directions 
(1,0),  (0,1)  (1,1),  each  occurring  with  multiplicity  n.  We  determine  all 
critical  points  of  the  polynomials 


Vx) 


Vj)e 


ijx 

, 


n  E  z+  . 
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SIGNIFICANCE  AND  EXPLANATION 


\ 

\ 

This  is  a  further  report  in  a  series  devoted  to  the  study  of  box 
splines.  Box  splines  have  been  Introduced  in  MFC  TSR  #2320  and,  provide  a 
natural  generalization  of  univariate  cardinal  splines,  i.e.,  splines  with 
uniform  knot  sequence. 

The  process  of  univariate  spline  interpolation  becomes  particularly 
simple  in  the  cardinal  case,  and  this  report  considers  the  corresponding 
bivariate  process  of  interpolation  at  the  integer  points  in  the  plane  to 
given  function  by  a  linear  combination  of  integer  translates  of  a  box 

spline.  In  particular,  the  report  shows  that  this  process  is  well  posed, 

C 

i.e.,  any  bounded  continuous  function  X  has  exactly  one  such  bounded 
interpolant  If.  The  argument  uses  the  Fourier  transform  to  identify  a 
certain  trigonometric  polynomial  (in  two  variables)  whose  nonvanishing  is 
equivalent  to  the  asserted  well-posedness.  The  minimum  value  of  this 
polynomial  yields  a  bound  on  the  norm  of  the  resulting  interpolation 
projector  I., 
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In  a  series  of  beautiful  papers,  I.  J.  Schoenberg  developed  the 
theory  of  univariate  cardinal  splines  [6-8],  A  basic  result  is  the 
positivity  of  the  Euler-Frobenius  polynomials  which  implies  the  well 
posedneas  of  cardinal  Interpolation. 

Theorem  1  [6].  Let  M_  denote  the  univariate  cardinal  B-spiine  with 


support  centered  at  0.  The  Euler-Frobenius  polynomials 


x  •  it  mod  2n  Z  (x  •  0  mod  2*  Z  ). 

In  this  note  we  obtain  the  bivariate  analogue  of  this  result  for 

box-splinea.  For  a  net  of  vectors  3  -  with  e  the 

box-spline  M-  is  the  functional  on  CgOt®)  defined  by  [ 1 1 


V  -  1,-1  l,n  «J,  W"' 


l  2»2j  1 


Equivalently,  M-  can  be  defined  by  its  Fourier  transform 
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where  S(z)  :■  (2/z)  sln(z/2).  The  latter  definition  stresses  the 
similarity  to  the  univariate  case.  We  define  the  multivariate  Euler- 


Frobenlus  polynomials  by 

(3)  P~(x)  l  M„(j)eijx. 

-  j^.m  a 


In  the  bivariate  case  (m  -  2)  we  proved  [3]  the  following  conjecture. 

The  polynomials  P_  are  strictly  positive  Iff  the  box-splines  M_(»  -  j), 

m 

j  e  Z  ,  are  linearly  independent. 

If  valid  in  general  (m  >  2)  the  conjecture  would  Imply  that  cardinal 
Interpolation  is  well  posed  if  the  obvious  necessary  condition  of  linear 
independence  is  satisfied.  For  two  variables  it  was  shown  in  [2]  that  the 
box-splines  are  linearly  Independent  only  on  the  "standard”  three 
direction  mesh,  up  to  symmetry  the  vectors  In  a  have  to  be  chosen  from 
the  set  {(1,0),  (0,1),  (1,1)}.  While  the  corresponding  grid  is  very 
regular,  the  analysis  of  the  interpolation  problem  is  complicated.  Our 
results  {3,4]  are  not  as  complete  as  in  I.  J.  Schoenberg's  univariate 
theory.  E.g.  we  were  not  able  to  determine  the  location  of  the  minimum 
for  P_.  which  in  general  depends  on  a.  We  conjectured  that  in  the 

symmetric  case,  when  each  of  the  three  vectors  in  a  occurs  with 
multiplicity  n,  the  polynomial  Pn  «  P~  attains  its  minimum  at  the  point 

(— ,y).  1°  this  note  we  prove  this  conjecture  and  determine  all  critical 

points  of  PQ. 


Theorem  2.  The  polynomials  P 


attain  their  minima  at 


(5“. j“)  ®od  2*  Z“ ,  their  maxima  at  the  points  2w  Z“  and  have  saddle 
2  2 

olnts  at  *  Z  \2a  Z  .  These  are  the  only  critical  points  of  P  . 


Figure  1  below  shows  the  gradient  field  of  Pj  on  [*/2,3*/2]  *  t~a,a] 
which  Illustrates  the  general  situation. 
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<  Figure  1  > 


The  proof  of  Theorem  2  relies  heavily  on  the  symmetries  of  PR.  Let 

A  denote  the  group  of  12  linear  transformations  which  leave  the  mesh 
generated  by  the  three  directions  (1,0),  (0,1),  (1,1)  invariant*  This 
group  is  generated  by  the  matrices 


~l  0 


0  -1 


1  0 


0  -1  j  ,  (-1  0  J  ,  U  -1 


-l  i 


0  1 


which  correspond  to  reflection  at  the  origin  and  permutation  of  the 
directions.  The  symmetric  box-spline  Mr  is  invariant  under  composition 

with  A,  i.e. 


M  (Ax)  -  M  (x),  A  e  A. 
n  n 


Therefore,  the  corresponding  Euler-Frobenius  polynomials  satisfy 


P  (A*x  +  2xj)  -  P  (x),  A  e  A,  j  e  Z  , 
n  n 


where  A*  denotes  the  transpose  of  A.  These  relations  give  much 
information  about  the  structure  of  P^.  Denote  by  Vf(u,v)  :■ 

(D^f (u,v) ,  Dvf(u,v))  the  gradient  of  a  function  f.  Differentiating  the 

Identity  (6)  we  obtain 


(VP  (A*x  +  2*j))  A*  =  Vp  (x),  A  e  A,  j  e  Z  . 


5 


Let  I  denote  the  unit  matrix.  Identity  (7)  Implies  In  particular  that 


(8)  Vp  (x)  e  ker(I  -  A)  If  (I  -  A*)x  -  2*j. 

n 

-1  0  f-1  1  ' 

For  A  •  ,  e  A  tt  follows  from  (8)  that  VP 

,  0  -1  J  (-1  0  J  n 

2  2ir  2ir  2 

vanishes  at  the  points  n  Z  and  +  2*  Z  respectively.  For 

1  -1  f-1  0  0  1  0  -1  f-1  1  1  0 

A  -  ,  ,  ,  ,  ,  e  A,  the 

,  0  -1  1-1  1  1  0  l-l  0  l  0  1  1  -1 

matrices  (I  -  A)  have  rank  one  and  it  follows  from  (8)  that  for  k  e  Z  , 

(0,l)7Pn(x)  -  0  If  (l,2)x  -  2*k, 


(l,0)VPn(x)  -  0  if  (2,l)x  -  2 irk, 


(9) 


(l,-l)VPn(x)  -  0  if  (l,-l)x  -  2*k, 
(l,l)VPn(x)  -  0  if  (l,l)x  -  2*k, 


(2,-l)VPn(x)  -  0  if  (1 ,0)x  -  2nk, 


(l,-2)VPn<x)  -  0  if  (0, 1  )x  -  2 irk. 


6 


The  remaining  4  matrices  in  A  give  no  further  Information. 

Let  ft  denote  the  (closed)  triangle  with  vertices  (0,0),  (*,0), 

(t”*t~)*  The  set 


ft*  -  U  Afl. 

Ac  A 

which  is  the  convex  hull  of  the  six  points  ±(-|^-,-|^),  ±(y^ ,  y^-),  ±(y~»  y-)» 

is  a  fundamental  domain,  i.e.  its  translates  form  an  essentially  disjoint 
2 

partition  of  3R  .  Therefore,  to  complete  the  proof  of  Theorem  2,  it  is 
sufficient  to  show  that 


% 


(10)  VPn(x)  *  0  for  x  e  ft  \  {(0,0),  (*,0),  (-f^*^)  1 
and  that 

(11)  Pn(T'*TL)  <  Pn(n’0)  < 

To  this  end  we  prove  the  following  estimates: 

(l)  DuPn(u,v)  /  (2u  +  v)  <  0  for  (u,v)  e  ft^  :=  {(u,v): 

0<^v£u,2u  +  v£  u  >  0}, 

(ii)  D^P^u.v)  /  (2x  -  2u  -  v)  <  0  for  (u,v)  e  ft^  :*  {(u,v): 

_<  2u  +  v  <  2x,  0  v,  u  +  2v  _<  2x}, 

(iii)  DvPn(u,v)  /  v  <  0  for  (u,v)  e  ft^  ;=  f(u,v):  0  <  v  =»  2x  -  2u  <  jj. 


f 


8 


2 

where  A  :■  2"  Z  .  For  (u,v)  e  U  and  large  n,  the  terms  with  | k. j  +  |  i 
small  dominate  in  the  expression  for  P^.  This  fact  is  crucial  for  the 

subsequent  estimates. 

Proof  of  (i).  We  write 


(13) 


DuPn(u’v)  ’  n  £  \  A  l 

(k,i)cA 


with 


afc  l  S(u  +  k)n_1  S(v  +  £)n  S(u  +  v  +  k  +  £)n" 

(14) 

^  S' (u  +  k)  S(u  +  v  +  k  +  £)  +  S(u  +  k)  S’(u  •>  .  +  k  +  £). 


Using  the  inequalities 


(15)  S(w)  ,  S’(w)  <  min  (1,2/  w  ), 


-  jj  w  <  s'(w)  <  -  w,  0<w<x, 


for  (u,v)  e  ilj ,  we  obtain  the  estimates 


bo,o  i  ‘  72  (2u 


+  v). 


t 


9 


b0,0  1  "  T6  u  S°°  ~  T6  (u  +  v)  s(tL)  i  ’  h (2u  +  v)* 

(16) 

j  bkt  ii  ,  4  »  8tt  r  8in(u/2)  ,  ain(  (irfv)/2 )  \ 

|bg  —  u+k  u+v+k+i  '  2u  +  v  2u  +  v  ' 


1 

6n 

prj 

u+v+k+i 

"  • 

For  (u,v)  e  and  (k,i)  *  (0,0),  we  have 


u+k  u+v+k+i  '  |v+l|  — 


-2 


Combining  this  Inequality  with  (16),  we  see  from  the  definition  of 
and  S  that 


,  Vn(u»v) 

'  nao,obo,o 


<  1 
A\(0 ,0 ) 


The  last  right-hand-side  is  less  than  1  for  n  ^  5.  Therefore, 
inequality  (i)  follows  from  the  second  inequality  in  (16)  and  the  fact 
that  ag,0  is  positive  on  Jlj. 


Proof  of 


into  the 


The  sets 


Therefore 


(17) 


where  (c. 


b 


m 


10 


(ii).  In  expression  (13)  for  Dupn  we  split  the  index  set  A 
three  parts 


Aq  :«  {(k,A):  2k  +  A  +  2*  -  o}, 

A+  :«  {(k,  A) :  ±(2k  +  A  +  2n)  >  o}. 


h+  and  A_  are  related  by  the  bijective  mapping 


(k.  A)  e  Af  <— >  (k' ,  A' )  *>  (-k  -  A  -  2ir,A)  e  A_. 


,  we  can  write  DP  in  the  form 
’  u  n 


Dupo(u..)  -  n  I  ^  ^  +  n  i  V.A.l 

\)  A+ 


f.  (14)) 


a-k-A-2Tr,  A 

,  A  °k,  A  afc  ^  -k-A-2*,  A 

=*  (S'(u+k)S(irfv+k+A)  +  S(u+k)S’ (irKr+k+A)  1 


+  5n-1  [S'(u-k-A-2ir)S(u+v-k-2it)  +  S(u-k-A-2n)S’(u+v-k-2w)] 


with 


11 


u  +  k  u+v»k»A 
u-k-A-2*  u+v-k-2* 


Observe  that  for  (u,v)  e  0^  and  (k,A)  e  Af 


(18) 


0  <  C  and 


2w-v-2u  2k+A+2* 

u-k-A-2*  ufv-k-2it  * 


Since  the  numerator  in  1  -  C  is  positive,  letting  A^:*  {(k,A)  e  Af: 
(k+A+ir)(k-Hr)  >  0},  we  have 


(19) 


0  <  C  <  1,  (k, A)  e  A*. 

0  <  i/C  <  1,  (k.  A)  e  A+\  A*. 


Using  the  identity 


(20)  s(p)s'(q)  ±  S'(p)S(q)  -  sin  8in  £  8in 

Pq  p  q 


we  can  simplify  the  above  expressions  for  b^  ^  and 


b.  .  and  obtain 

It  |  w 


(21) 


2  sin(u4'v/2-x) 

V, A  “  (u+k)(ufv+k+AT 


4(2ufWit)  .  .  8ln  u£  flln  wfvfktA  #  (kfA)  e  A 
(u+k)Z(ufv+k+A)Z 


12 


;  2(-l)*'  sln(trt-v/2)/,  .  4(-l ) 1  sin(u/2)sin((u+v)/2) 

<22)  V*  -  («k)(^«rU  (u*k)2(irt-W*Hl>2 


x  [(2irHv+2k+A)+Cn+1(2vrt-v-2k-A-4x)],  (k,i)  e  A*. 


In  the  terra  In  square  brackets  we  add  and  subtract  (2u  +  v  -  2k  -  i  - 
4x).  Then  a  direct  computation  using  (18)  yields 


(23) 


[...]  -  (2u  +  v  -  2*)(2  + 


(2k+fc+-2tt)(2u+v-2k-&-4tt) 
( u-k- 1-2 *)( u+v-k-2  X ) 


Analogous  to  case  (i)  we  show  that  a^  ^b^  q  is  the  dominant  term 
for  the  right  hand  side  of  (17).  Indeed, 

<24>  Vo  ± --6  lESSr  f"  "^5- 

as  one  checks  numerically  for  n  ■  5,  and  therefore  has  it  for  n  5, 
since  bQ  Q  decreases  as  n  increases  as  we  see  from  (19),  (22)  and 
(23).  For  (u,v)  e  we  have  x/3  £  u,  u+v  £  4ir/3  and  we  obtain  from 

(19 )— (24 )  the  estimates 


1  <  i_  |_u_| 

1  u+v  1 

1-  1  -  .6  lu+kl 

I  u+v+k+il  ’ 

0,0 


0t,*)  e  AQ, 


(25) 


13 


k. Ai  <  3(iH-l)(2k*&f2*)  u  i__u+v 
-  —  .6  o+k  u+v+k- 


o+k|  lu+v+k+il’  (k,A)  6  A*' 


For  (k,i)  e  A+\  A^  we  estimate  C  “b^  4 


In  a  similar  way  and  obtain 


k.li  „  n  3(n+l)(2kfl+2it)  u  i  i  u+v  i  /u  •>  r  *  v  » 

r-*i  i  c  re  ^+k|  |G+7+k+I|‘  (k**>  e  V  V 


For  (k,i)  "  (0,-2*)  we  obtain  the  sharper  estimate 


r^l  <  -‘  .  •>*. 


numerically  for  n  ■  5,  hence  valid  for  n  5  since  ^  increases 


with  n. 


Similarly  as  for  case  (i),  it  follows  from  (17),  (25)-(28),  the 
definition  of  C,  and  the  Inequality 


(v/(2*  -  v))n((u  +  v)/(2*  -  u  -  v))n  <  1,  (u,v)  e 


* V  *  » *  ,  1 
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D  P  (u,v) 


nao,obo,o 


-  1  <  .6  + 


2  r  i  it  in  i  5*/6  in  i  4*/3  m  . 

.6  I** I  15776+11  |4it/3+k+i| 


A0\(0,-2it) 


3(trf  1 ) 


rP.i  i*r ifer «**'> * 


A*\(0,0) 


fe  ^  ifeiin  ifer  ifer  (2k+«->' 


The  right-hand  side  is  less  than  1  for  n>  5  and  the  inequality  (ii) 
follows  from  (24)  and  the  fact  that  ap  q  is  positive. 


Proof  of  (ill).  We  have 


W*  -  v/2‘ v)  - n  {  K,iK,i 


a'  e  S(ir  -  v/2  +  k)n_1S(v  +  i)n_1S(ir  +  v/2  +  k  +  i)n  1 


b£  £  S(u  -  v/2  +  k)S'(v  +  *)S(*  +  v/2  +  k  +  i)  + 


S(n  -  v/2  +  k)S(v  +  i)S'(*  +  v/2  +  k  +  i). 


Note  that  aQ  Q  -  *_2*,0'  lt  can  be  verlfied  n«nerically  that 

(30)  C  s-  -  sup  (b*  Q  +  b’  0  )/v  >  .1. 

0<v<n/3  '  f 

To  estimate  the  remaining  terms  in  (29)  we  observe  from  the  definition  of 
S  and  (15)  that 

1  2v2  min  t1,  | n-v/2+k|  |v+l|  I ir+v/2+k+ij' 

Therefore, 


The  right-hand  side  is  less  than  1  for  n  _>  2  which,  together  with 
(30),  implies  the  inequality  (iii). 
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